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ABSTRACT 



Context. We present evidence that parsec-scale jets in BL Lac objects may be significantly distinct in kinematics from their counter- 
parts in quasars. 

Aims. A detailed study of the jet components' motion reveals that the standard AGN paradigm of apparent superluminal motion does 
not always describe the kinematics in BL Lac objects. We study 0735+178 here to augment and improve the understanding of the 
✓ ■ peculiar motions in the jets of BL Lac objects as a class. 

Methods. We analyzed 15 GHz VLB A (Very Long Baseline Array) observations (2cm/MOJAVE survey) performed at 23 epochs 
between 1995.27 and 2008.91. We then compared the jet kinematics with the optical and radio light curves available for this BL Lac 
object and point out some striking correlations between the properties of the radio knots and the features in the light curves. 
Results. We found a drastic structural mode change in the VLBI jet of 0735+178, between 2000.4 and 2001.8 when its twice sharply 
bent trajectory turned into a linear shape. We further found that this jet had undergone a similar transition sometime between December 
1981 and June 1983. A mode change, occurring in the reverse direction (between mid-1992 and mid-1995) has already been reported 
in the literature. These structural mode changes are found to be reflected in changed kinematical behavior of the nuclear jet, mani- 
fested as an apparent superluminal motion and stationarity of the radio knots. In addition, we found the individual mode changes to 
correlate in time with the maxima in the optical light curve. 

Conclusions. We found for this BL Lac object a drastic change in the kinematics of the nuclear jet, i.e, transition from 'typical super- 
luminal' to an unusual 'stationary' state. Interestingly, we found that this change is accompanied with a mode change in the nuclear 
jet's morphology. The long sequences of VLBA images reported here for this BL Lac object do not support the commonly assumed 
■ connection between radio flux-density outbursts and the ejection of new VLBI components. 

£C) • Key words. BL Lacertae objects: individual: PKS 0735+178 - Radio continuum: galaxies - Techniques: interferometric 
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f^) , 1. Introduction tudes are typical for radio selected BL Lacs (Gopal-Krishna et 

t-H . al. 2003). Gupta et al. (2008) claimed detection of optical vari- 

\ The strong radio source PKS 0735+178 is an optically bright ability with an amplitude of 4.9% on one night, but since that 

• i-H ! classical BL Lac object, highly variable at radio/optical wave- light curve remains steady, any such variability must be of the 

^ ■ lengths (e.g., Goyal et al. 2009 and references therein). Based flicker type and not on hour-like time scale. 

H | on an absorption feature in the optical spectrum, a lower redshift In the radio band, an unusual property highlighted for this BL 

. . . , limit of 0.424 was established (Abraham, et al. 1991; Stickel et Lac is its remarkably flat radio spectrum, termed as "cosmic con- 

al. 1993; Scarpa et al. 2000). This is consistent with the non de- spiracy" (Cotton et al. 1980). 
tection of the optical host in the HST snapshot imaging survey 
by Sbarufatti et al. (2005). A companion (possible absorber) has 

been detected at a projected separation of 22 kpc (Pursimo et al. 1.1. Radio morphology: kpc- to pc-scale 
1999). Extensive monitoring of this blazar in the radio and op- 
tical bands revealed variability on different timescales. Periodic PKS 0735+178 has been observed with radio interferometers at 
features of 13.8-14.2 years were found by applying the Jurkevich various frequencies. The 1.64 GHz VLA (Very Large Array) ob- 
method to about a century long optical monitoring data (Fan et servations by Murphy et al. (1993) revealed a rather compact 
al. 1997; Qian & Tao 2004). Ciprini et al. (2007) reported possi- structure, slightly elongated at a position angle of ~ 160 deg. The 
ble detection of several shorter time scales. In their search for 5 GHz MERLIN (Multi-Element Radio Linked Interferometer 
intra-night optical variability on 17 nights between 1998 and Network) image taken in 1983 is dominated by a core and also 
2008, Goyal et al. (2009) did not detect on any occasion intra- shows a very weak component to the southeast (Baath & Zhang 
night variability with an amplitude >3%, though such ampli- 1991). The core-to-component flux density ratio is 200:1 and the 
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position angle of the component is ^160 deg - in agreement with 
the VLA observations by Murphy et al. (1993). These observa- 
tions established a large misalignment between the kpc-scale and 
pc-scale structures of this source by ~ 90-100 deg. 
PKS 0735+178 has been a favorite target of VLBI monitoring 
programs, and striking changes in its nuclear jet are revealed 
from the more than two decades of its VLBI monitoring. 
Baath & Zhang (1991) imaged the source at 5 GHz with global 
VLBI arrays at four epochs between 1979 and 1984. The large 
structural changes thus found were interpreted by them in terms 
of relativistic (y ~ 8) motion of the radio knots along a curly 
path with a steadily decreasing viewing angle. Particularly in- 
structive are their VLBI images of the inner 5 mas region, taken 
at the epochs 1981.9 and 1983.5, which (despite modest resolu- 
tion) reveal a clear changeover from a kinky staircase profile in 
1981.9 to a straight mode in 1983.5 (see below). This morpho- 
logical transition, first noted in the present work, is discussed in 
Sect. [3] 

Gabuzda et al. (1994) observed the source with the VLBA at the 
epochs 1987.41 (5 GHz), 1990.47 (8.4 GHz), and 1992.44 (5 
GHz). In all these observations, the source revealed a straight jet 
extending to the northeast. Gomez et al. (2001) present VLBA 
observations at 8.4 and 43 GHz between August 1996 and May 
1998. Between March 1996 and May 2000 Agudo et al. (2006) 
observed the source using the VLBA at frequencies between 5 
and 43 GHz. The observations by Gomez et al. and Agudo et 
al. show a twisted jet with two sharp apparent bends of 90deg 
within two milli-arcseconds of the core. Below we call this tra- 
jectory staircase mode. Between mid-1992 and mid 1995, the 
jet trajectory changed again drastically from straight to stair- 
case (Gomez et al. 2001). Unfortunately, no VLBI observations 
documenting this second transition are available. 

From early VLBI observations of this BL Lac object Baath 
& Zhang (1991) estimated bulk Lorentz factors of up to ~ 8/i _1 c 
for the knots in the VLBI jet, moving from the core towards a 
stationary component located ~ 4 millisecond northeast of the 
core. A subsequent paper by Gabuzda et al. (1994) reports an 
apparent superluminal motion of >7.4 hr l c of one knot. After 
mid- 1995 the components show a tendency to cluster near sev- 
eral jet positions, remaining stationary in time (Sects. I3land l4.2l i. 

1.2. The present work 

The extensive VLBA coverage of this blazar at 15 GHz during 
1995-2009 under the MO J AVE program allowed us to undertake 
a detailed probe of structural transitions in the nuclear jet of 
this BL Lac object (Sect. |3). Thus, in addition to the last 
documented mode change in the pc-scale trajectory (Gomez et 
al. 2001), we found clear evidence for another mode change 
between 2000.4 and 2001.8. We find that from the staircase 
mode observed since mid 1995 the jet trajectory has reverted to 
a straight morphology during the time interval covered by our 
VLBA data analysis, so that the timing of this morphological 
change can be pinned down to an uncertainty of merely ~ 1.6 
years. In comparison, the uncertainty for the mode change 
inferred from the data of Gomez et al. was ~ 3 years. According 
to them, the morphological transitions are most likely caused 
by pressure gradients in the external medium through which the 
jet propagates. We discuss this hypothesis with regard to the 
mode transitions found by us. We further examine these mode 
transitions gleaned from the VLBI data, as mentioned above, 
with regard to the observed trends in the radio and optical 
light-curves of this intensively monitored blazar. 



This paper is organized in the following way: we first present 
the results of our analysis of the archival VLBA data (MOJAVE) 
and point out the newly noticed correlations between the mor- 
phological and flux-density variations of this source. We then 
discuss our findings in the context of previously reported infor- 
mation about this blazar. In particular, we explore the connection 
between the peculiar mode changes and kinematics of its nu- 
clear jet. Apparent transverse velocities (/3 app ) in this paper were 
calculated with the one-year WMAP (Wilkinson Microwave 
Anisotropy Probe) data (Spergel et al., 2003); the values adopted 
are h = 0.71+°$; Q. m h 2 = 0.135!°°°*; Q tot = 1.02 + 0.2. The 
one-year and three-year WMAP parameters are less than 0.1c 
for n = 0.1 masyr -1 out to z = 4; this is small with respect to 
the formal uncertainties in the apparent velocities. 

2. Data analysis 

We analyzed the VLBA observations for 23 epochs, taken at 15 
GHz under the 2cm/MOJAVE survey between 1995 and 2009. 
The data were model-fitted with the difinap model-fitting rou- 
tine. To each epoch's VLBA image we attempted to fit the min- 
imum number of components needed to describe the structure 
adequately to facilitate a least ambiguous component identifica- 
tion from epoch to epoch. For each image the model-fitting pro- 
cess started with a point-like source model. In order to keep the 
comparison of the images for different epochs simple and credi- 
ble, we assumed the components to be circular; highly elliptical 
components are not likely anyway. 





Fig. 1. a) Core separation of each VLBI component as function of time for both time intervals described in the text. Grey symbols 
represent weak components that could not be identified through the epochs, b) Position angle of each component as function of 
time, for both time intervals, c) The jet ridge-line for all the epochs observed, d) Three representative jet ridge-lines selected from 
the two time intervals, for comparison, e) Flux-density evolution of the individual components and the core (upper data points), f) 
Radio flux-density light-curves from the Michigan monitoring program (UMRAO). Indicated are the two transition phases where 
the mode changes (i.e., morphological transitions) occurred (Sect. 3). 
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1980 1985 1990 1995 2000 2005 201 

Epoch of measurement 

Fig. 2. The diagram, reproduced from Goyal et al. (2009), shows 
the multiple light-curves of PKS 0735+178. Also displayed 
are the R-band light curve (a), the B-Band light curve (b), the 
UMRAO light curves at 15 and 5 GHz (c) & (d) and the run of 
the radio spectral index (e). The plot (f) shows the variation of 
the percentage (linear) polarization at 15 GHz. The two pairs of 
vertical lines bracket the structural transition phases discussed in 
the text. 



our observations. In Fig. [6] the jet shows the twice bent struc- 
ture - the staircase mode, whereas in Fig. [7] the jet structure is 
seen to have become straight. Figure [1] a) shows the variation of 
core separations between the epochs 1995.27 and 2008.91 for 
the different components. Striking variations in the jet kinemat- 
ics are seen during the period 1995 - 2000.4. Four jet compo- 
nents clearly show outward directed motion with apparent super- 
luminal speeds, while one component shows an inward motion 
(directed towards the core) also with an apparent superluminal 
speed. The apparent speeds derived from linear regressions to 
the component positions are given in Table Q] From 2001.8 on, 
only four components can be identified in the VLBI jet, and the 
jet components are found to be stationary (Figs.|6]and|7]). 
Figure fl~|b) presents the position angle (relative to the core) at 
different epochs for each component. Again we find a significant 
change to have occured between 2000.4 and 2001.8. Between 
1995 and 2000 the position angle distribution is quite broad with 
values ranging from ~10°-100°, which amounts to a wide spread 
of over ~90°. After 2000.4, the position angle distribution is seen 
to become much narrower, with a total spread of just about 20° 
around a mean value of ~65°. 

In Fig. Q] c) the relation between the core separation and 
position angle is plotted for each of the 23 epochs. This plot 
more clearly depicts the two contrasting trends present in the 
VLBI images. Thus, a different pattern is observed for the 
epochs before 2000.4 and after 2000.4. Prior to 2000.4 all the 
jet ridge lines show a similar form, with a clear minimum at 
a core separation of ~1.5mas, where the position angle drops 
to ~ 20°. After 2000.4, the form of the jet ridge lines becomes 
almost linear. 

Thus, both the jet morphology and kinematics are found to 
change remarkably between 1995 and 2009, as summarized 
below: 



1. In the staircase jet mode five radio knots/components move 
with apparent superluminal speeds either away from the 
core (four components), or towards it (one component). The 
position angle distribution is broad (~ 90°). The jet ridge 
lines for all epochs before 2000.4 are similar, with a clear 
minimum in position angle. 



2.1. Estimation of the uncertainties 

Although the use of VLBA data generally ensures a high qual- 
ity, considerable variation can occur between the epochs in both 
the amount of data and its quality. The uncertainties were de- 
termined individually for each epoch by comparing the model- 
fit parameters for the sets with different numbers of components 
(+1). For the uncertainties of the flux-density we assumed a min- 
imum error of 5%. The positional uncertainty was determined 
by averaging l/5th of the beam size and the value determined by 
comparing several model-fit results. For the size we adopted an 
uncertainty of 10%. The model fit parameters and their uncer- 
tainties are given in Tables [3] and [4] We also list there the epoch 
of observation, the flux-density of each fitted component and its 
separation from the core, the corresponding position angle, and 
the component size. 

3. Results 

Figures|6]and|7]show the images derived by us, with the model- 
fit components superimposed. It is evident that the jet's morphol- 
ogy undergoes a drastic change during the time span covered by 



2. In the straight jet mode, separations of all four VLBI 
components from the core remain essentially unchanged; 
their apparent speeds are sub-luminal. The position angle 
distribution is narrow (~ 20°). 

In Fig.[T]d) we present the same information as in Fig.[T]c), but 
for selected epochs only: the two epochs before the mode change 
and the first three epochs following the mode change (2001.81, 
2001.98 and 2002.90). Here the different jet ridge lines for the 
two modes are seen more clearly. FigureQ]e) displays the evolu- 
tion of flux-densities of the individual components and the core; 
again significant changes are observed. 

3. 1 . Morphological transitions and flux-density variability 

The BL Lac 0735+178 is known for its pronounced flux variabil- 
ity in both optical and radio bands (Ciprini et al. 2007; Goyal 
et al. 2009). The origin of this variability, however, is still de- 
bated. To check for a possible correlation between the opti- 
cal/radio flux-density variability and the morphological transi- 
tions witnessed in the VLBI observations, we show the light 
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Table 1. Proper motion and apparent speeds of the individual jet 
components in 0735+178 (Sect. 3). 



mode 


C omponcnt 


fi [mas/year] 


B Tel 

Hupp L U J 


staircase 


C 1 


U.Uj ± U.U1 


n 11 _i_ n ia 
U. / / ± U.zo 




c2 


0.46 ±0.16 


11.87 ±4.16 




c3 


0.22 ± 0.02 


5.68 ± 0.52 




c4 


0.22 ± 0.04 


5.68 ± 1.03 




c5 


-0.13 ± 0.03 


-3.36 ± 0.77 


straight 


ca 


0.01 ±0.01 


0.26± 0.26 




cb 


0.02 ± 0.02 


0.52 ± 0.52 




cc 


0.01 ± 0.02 


0.26 ± 0.52 




cd 


0.01 ± 0.02 


0.26 ± 0.52 



curves (Figs. |T]f, |2j) with the epochs of the observed mode tran- 
sitions indicated by pairs of grey and black lines, respectively. 

3.1.1. Radio flux variability 

In Fig. Q]f) the intergrated flux-densities at 4.8, 8.0, and 14.5 
GHz are plotted, monitored under the University of Michigan 
programme (UMRAO). The grey stripes denote the two most re- 
cent mode transition events. Based on the VLBI images reported 
here and in the literature (see above), one can deduce that the first 
of these mode transitions (straight to staircase) occured some- 
time between mid-1992 and mid-1995, and hence its time-scale 
must be less than three years. Likewise, the latest mode transi- 
tion {staircase to straight) we report in this paper took less than 
two years. The time interval between the two transitions is ~7.5 
years. 

In Fig. [2] we show a plot adopted from Goyal et al. (2009) 
and indicate in it all three known mode transitions within pairs of 
vertical black lines. The 2001 transition occured shortly before 
a clear radio flare. The 1994 transition occured shortly before 
a clear modest radio flare. The 1983 mode transition cannot be 
related to a significant radio flare. While the corelation between 
the times of the transititons and the optical flaring is convincing, 
the case is not so clear with regard to a possible correlation 
with radio flares. It seems that for the last two transitions the 
radio flares are delayed with regard to the optical flare, which 
could be explained by propagation effects within the source. 
Unfortunately only one VLBI observation - the observation by 
Gabuzda et al. (1994) from 1990.47 - is available for the time 
of the most prominent radio flaring around 1990. This would 
help to understand the relation between the radio flaring and 
the morphological changes. A simultaneous monitoring of the 
radio and optical flux in combination with multi radio frequency 
VLBI observation is required for a better understanding of the 
emission processes. 



3.1.2. Correlation with optical flux variability 

In Fig. [2fa, b) we reproduce the optical light curves from Goyal 
et al. (2009), on which we have marked the three known mode 
transitions mentioned above. Both staircase to straight mode 
changes clearly take place near the brightness maxima. The op- 
tical flares coincident with these two transitions are more promi- 
nent compared to the middle flare coincident with the opposite 
mode transition. We thus find that a strong correlation exists be- 
tween the structural change in the parsec-scale jet and the optical 
light curve. To our knowledge, such a correlation has not been 



pointed out earlier for any jet, and this underscores the impor- 
tance of long-term, coordinated VLBI and optical monitoring of 
blazars. 

The origin of this conspicuous correlation in this source needs 
to be investigated in greater detail. This correlation could be ex- 
plained by a small change in the direction to the line of sight. 
This would change the amplification of intrisinc jet bending and 
also allow for a more direct view into the inner regions of the 
AGN. This "geometric" explanation should be explored further. 
While more subtle interpretations are possible, a simple possi- 
bility to understand this correlated radio-optical behavior is that 
the optical flare manifests an impulse by which a large amount of 
kinetic power (momentum flux) is injected by the central engine 
into the nuclear jet, as a result of which the jet gets stretched 
into a linear shape. This explanation is related to the question 
of higher stability of powerful jets and the collimation process 
of jets in general (e.g., Nakamura & Meier 2009). The observed 
time lag between the optical flare and the disappearance of the 
jet's kink could in fact be quite small (< 1 yr), since the jet flow 
in blazars is directed close to the line-of-sight, which is also the 
direction of propagation of the optical signal recorded. An al- 
ternative scenario is presented in Sect. 14.41 Either of these two 
possible models will also be required to explain the disappear- 
ance of superluminal motion following the mode transition to a 
straight jet. 

The proposed energetic impulse that straightens the kinky jet 
could originate from the core. Assuming that 0735+178 is a heli- 
cal jet, whose components are therefore seen at different viewing 
angles at the times of different mode transitions, then this could 
explain the connection to the radio flares as well. The radio flare 
would then not be caused by the ejection of a new component - 
which is usually assumed - but by the Doppler factor changing 
with the viewing angle. 

4. Discussion 

4. 1 . Morphological mode transitions 

PKS 0735 + 178 showed a kinky jet trajectory in the VLBI 
observations in December 1981 and a straight jet trajectory 
in June 1983. Therefter, sometime between mid-1992 and 
mid-1995, the morphology changed back to kinky with two 
sharp bends of ~ 90° (staircase mode), within 2 msec of the 
core. Between 2000.4 and 2001.8 the jet trajectory changed 
once again to the straight mode, as found in the present work. 
Thus, three clear mode changes of the pc-scale trajectory of this 
jet were observed. The time gaps between these transitions are 
about 11.3 years and 7.2 years. A simple scenario for the mode 
transition from kinky to straight jet was outlined in the previous 
section, taking a clue from the observed coincidences of the 
morphological transition events (kinky to straight jet) with the 
optical flares of this blazar. Below we present a kinematic model 
for the latest mode transition revealed by our analysis of the 
VLBI observations taken at 23 epochs between 1995 to 2009 
(Sect. 1441). 



4.2. Jet kinematics 

During the staircase mode, outward and inward apparent 
superluminal motion is observed. During the straight mode, 
the separations of various knots/components from the core 
remained essentially unchanged. Whether these components are 
the same or different, cannot be decided with absolute certainty. 
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However, should the mode transition be caused by a changing 
viewing angle - a possibility argued below - the components 
could be the same. 

A similar phenomenon - rather stationary components relative 
to the core, but moving perpendicular to the jet ridge line - was 
studied by us in the two prominent BL Lac objects, 1803+784 
(Britzen et al., 2009) and 0716+714 (Britzen et al., 2009). 
Over a time of ~ 20 years we do find rather slow apparent 
motion or stationarity for nearly all of the components in the 
pc-scale jet of 1803+784. However, the components are seen 
to move significantly in position angle measured from the core. 
The entire ridge line of the jet thus evolves in the transverse 
direction, changing from a straight line to an almost sinusoidal 
form and back again to a straight line. These significant jet 
ridge-line changes occur on a timescale of ~ 8 years. A similar 
phenomenon was studied by us in the highly variable BL Lac 
0716+714. Here again we find rather stationary components 
that do, however, show motion in position angle and, once 
again, the jet's ridge line evolves. We expect that the stationary 
components can be explained in terms of a geometric origin. In 
the past, various mechanisms were proposed to explain slow or 
subluminal apparent motions. Smaller apparent speeds can be 
obtained by a flux-weighted integration of the apparent motions 
over the entire jet width (Gopal-Krishna et al. 2004). For a brief 
review we refer to Britzen et al. (2008) and references therein. 
Different studies have reached opposite conclusions on whether 
BL Lac objects typically have slow or fast jet components 
(e.g., Gabuzda et al. 1994; Wehrle et al. 1992; Ghisellini et 
al. 1993; Jorstad et al. 2001). The CJF (Caltech-Jodrell Bank 
Flat-Spectrum) survey of 293 AGN (Britzen et al. 2008) shows 
some tendency for smaller apparent velocities of BL Lac objects 
compared to quasars. However, the small size of the BL Lac 
sample handicaps the comparison. 



4.3. A "hybrid" between quasars and BL Lac objects? 

As noted in the previous subsection, there has been an ongoing 
debate in recent years whether BL Lac objects and Quasars re- 
veal similar kinematic properties or whether the apparent speeds 
differ significantly. According to the general picture of the AGN 
paradigm, fast apparent speeds with values up to ~30c can be 
observed (e.g., Piner et al. 2006) and in almost all cases the mo- 
tion is directed away from the assumed stationary core. This is 
in contrast to what we find in a growing number of BL Lac ob- 
jects. Based on our analysis of 0735 + 178 presented here and the 
detailed studies of 1803+784 and 0716+714, it is possible that 
the kinematic properties are significanlty different. 0735+178 is 
a "hybrid" with regard to the mode changes and the accompany- 
ing kinematics - it shows 'quasar'-like fast apparent motions at 
certain periods of time and then again 'slow' motion (as we find 
in BL Lac objects) at other times. Thus this source switches be- 
tween the different states that we tend to associate with quasars 
or BL Lac objects. Due to this atypical behavior, 0735+178 can 
help us to understand the difference between these two types of 
AGN. It is often assumed that the basic difference between these 
two classes is based on a different line of sight to the object, 
with the BL Lacs being seen under the smallest angle to the line 
of sight (Unified Schemes). Within this scheme, it is then quite 
conceivable that the mode changes are caused by small changes 
in the angle to the line of sight, e.g., the jet bends slightly to- 
wards or away from us. Undoubtedly, 0735+178 is an extreme 
object since these changes of the angle can not be very consid- 
erable in the source frame. Thus, 0735+178 might already be 



da/at 




Fig. 3. This figure is reproduced from Gong (2008) and shows 
the difference between ballistic and non-ballistic motion. The 
model of ballistic motion plus precession is represented by the 
solid spiral connecting the features a, b, and c. Denoted by the 
dashed ellipse through feature b is the non-ballistic model which 
precesses with a constant distance to the core. 

oriented in a critical angle towards us. This would explain why 
this is the only source so far where these mode fluctutations have 
been observed. 

4.4. Helix or pressure gradients in the external medium? 

Gong (2008) proposed a non-ballistic model (see Fig. [3}, in 
which a continuous jet produces discrete hotspots via the in- 
teraction of the relativistic outflow with the ambient medium. 
Thereby significant bulk energy of the jet is transferred to the 
swept up material. The precession of the jet axis results in the 
motion of hotspots at different directions along the precession 
cone, which resembles a hotspot moving along a ring-like tra- 
jectory. Significant bulk energy loss may occur repeatedly for 
a continuous jet with a fixed axis, leading to several distinct 
hotspots at different distances from the core. The precession of 
such a continuous jet results, on the other hand, in the motion of 
hotspot at different rings and with different precession phases (a 
detailed analysis is in preparation). 

Under this scenario, the rectangular structure is just a state 
of the motion of these hotspots in different rings, which can ap- 
pear and disappear depending on the precession phases of these 
hotspots and their projections on the sky plane. Thus, conceiv- 
ably, the non-ballistic model of Gong (2008) can provide a fairly 
simple explanation for the appearance of the staircase mode. 

Figure [4] shows the observation of the BL Lac object 
0735+178 at 8.4 GHZ by Gomez et al. (2001) as well as the 
curves fitted to the data. The inner components K8, K6 and K5 
- as defined in Gomez et al. (2001) - are fitted within the non- 
ballistic model (Gong 2008) by three solid curves (correspond- 
ing to ring 1 , ring 2 and ring 3, respectively) at different distances 
to the core. The evolution of the three curves is clockwise, all of 
them start in 1988.47 and end in 2002.97. 

Under the simplest treatment, the non-ballistic model just 
gives the same opening angle for the precession cone for the 
three rings, A\ = X% = A3. However, in this case the morphology 
of PKS 0735+178 cannot be fitted well. 

Thanks to the scenario of the pressure gradients in the ex- 
ternal medium through which the jet propagates (Gomez et al. 
2001), it is reasonable to assume that A\, A%, A3 can deviate from 
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Fig. 4. The observations of the BL Lac object 0735+178 at 8.4 
GHz from Gomez et al.(2001) as well as the curves fitted to the 
data. 

Table 2. The parameters applied to fit the data of 0735+178 in 
the non-ballistic model. 



n 


i 




A, 


X 2 


^3 


15.0 


9.8 


225.6 


3.5 


3.0 


4.6 


1i 


12 


13 


Ri 


R 2 


R 3 


165.2 


21.7 


-53.0 


14.0 


11.1 


4.5 



each other by a few degrees. Under these circumstances, the evo- 
lution of the staircase mode can be well fitted, as shown in Fig. [4] 
The input parameters to this model are given in Table [2] 

The precession speed, Cl, is in units of deg/yr, the distance 
of a hotspot to the core, Ri, is in units of mas, and all other pa- 
rameters are in units of degree, i is the inclination angle between 
the jet's rotation axis and the line of sight; £ is the coordinate 
transformation angle; A, is the opening angle of the precession 
cone; 77, is the initial phase. 

Although Fig.[5]shows the fit to the 8.4 GHz observations of 
Gomez et al. (2001), the theoretical trajectory (Fig. [5]h) can be 
extended to up to 2009.63. By doing this, the disappearance of 
the staircase observed at 15.335 GHz and 15.365 GHz (symbol- 
ized by the ellipses) at the end of 2001 can also be explained as 
shown in Fig. [3] This figure shows the evolution of the jet ridge 
line based on the observed morphology (filled filled squares and 
ellipses) and the model (open squares). For the time of 2009.637 
the model prediction is shown in the last panel of this figure. 
This figure distinctly shows the different morphological phases 
of the jet ridge line as they are observed and how well they can 
be reproduced within the model. 

Consequently, the staircase mode occurs when the hotspot 
precesses to certain phases so that the projected appearance is 
rectangular. Likewise, the staircase mode disappears when the 
hotspot precesses to other phases where the projection devi- 
ates from the rectangular image. Obviously, the precession of 
a hotspot varies the angle of the jet axis from the line of sight. 
This leads to a varying Doppler boosting, which explains the 
observed radio flux-density changes. The non-ballistic model, 
as outlined before, can explain the observed phenomena. When 



interpreted in the context of the ballistic model, the two sharp 
apparent bends of 90 deg should move outward as a whole rel- 
ative to the core. Since each component in the staircase mode 
corresponds to a knot ejected previously, it should not have re- 
mained stationary. The staircase mode may induce pressure gra- 
dients in the medium through which the jet propagates, trig- 
gered by its own interaction with the ambient medium. In this 
case, gradual changes in the position and curvature of the jet 
bends near K8 and K6 with time can be expected (Gomez et 
al. 2001). Plausibly, such a scenario could be responsible for 
the small angle variations of the jet. However, to explain the 
two 90 degree bends would require very special conditions in 
the ambient medium, making it hard to understand why such 
an ambient medium existed only during 1992-2000, but was 
absent before mid- 1992 and again after 2001. As pointed out 
by Gomez et al. (2001), the moving-component scenario sug- 
gests very non-ballistic trajectories for K5 and K6. But there are 
still unanswered questions concerning the origin of these non- 
ballistic trajectories. Are they caused by continuous or discret 
ejection events, and what role, if any, is played by the jet preces- 
sion? Further detailed investigations and modeling is required to 
answer these questions. 



4.5. Timescales: Flux-density variability and mode change 

The three drastic mode changes are clearly correlated with 
maxima in the optical light-curve. Thus we assume that the 
optical emission (most likely emergent from the core) causes 
the change in the trajectory. In two cases, a major long-term 
radio flux outburst occurs shortly after a mode change has 
appeared and this can be explained within the non-ballistic 
motion scenario of hotspots outlined in the previous subsection. 
In the optical emission, Fan et al. (1997) and Qian & Tao (2004) 
found periodic features of 13.8-14.2 years. This is roughly twice 
the time interval we find here between the two mode transitions. 
A physical connection between the morphological changes 
(mode changes) and the flux variability in the radio and the 
optical regime is thus likely. 

According to conventional wisdom, radio outbursts are ex- 
plained as correlated with the ejection of a new radio component 
from the core. However, in PKS 0735+178 we found a similar 
number of components over a time span of roughly 14 years and 
no new components were unequivocally detected during this 
period. In addition, we find a correlation with the morphological 
changes described as mode changes. Furthermore, we are able 
to explain the flux-density changes within the model by Gong et 
al. (2008) of non-ballistically moving hotspots. Thus, it can be 
inferred that in PKS 0735+178 ejection of VLBI components 
from the core is not required to explain the observed (quasi- 
periodic) radio flaring. 

We found similar results for the other well-known blazars, 
1803+784 and 0716+714 - the long-term radio variability can 
be correlated with morphological changes in the nuclear jet and 
is hence also likely to be of geometric origin. We thus find in 
three BL Lac objects supporting evidence for a correlation be- 
tween long-term radio variability and morphological changes 
that are likely to be of geometric origin. Additional contribution 
to morphological changes, particularly the more abrupt ones, 
could come from short bursts of power injected into the jet by 
the central engine, manifested by the optical flaring events. 
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Fig. 5. The observed morphology of 0735+178 is indicated by filled squares and ellipses, corresponding to the VLBI observations 
at 8.4 GHz (1990-1998) and at 15 GHz (2000-2007). The evolution of the trajectory of hotspots (hollowed square) from 1990 to 
2009 is reproduced in terms of the non-ballistic model. 



5. Conclusions 

- We report a third (and the most clear) event of structural 
mode transition in the nuclear jet of the BL Lac object PKS 
0735+178. In particular, we find that the twice sharply bent 
trajectory of the nuclear jet transformed into a straight jet 
sometime between 200.4 and 2001.8. 

- From a literature search we further found that a similar tran- 
sition from a kinky to a straight jet is present in the VLBI 
observations of this source performed at 5 GHz in December 
1981 and June 1983. 

- Contrasting kinematic properties are found associated with 
the two morphological modes of the nuclear jet: in the stair- 
case mode apparent superluminal motion of the knots in the 
jets is observed, either directed away from the core or to- 
wards it. In the straight mode, however, the knots remain 
fairly stationary relative to the core. Their position angle dis- 
tribution is broad in the former case (~ 90°) and narrow in 
the latter (~ 20°). 

- The time interval between the mode changes is ~ 11.3 years 
between the first two mode changes and ~ 7.2 years between 
the last two mode changes. 



- Both mode changes, from kinky to straight jet, are found to 
coincide with individual, well-defined peaks in the optical 
light-curve. This points to a causal connection between the 
optical flaring and the structural changes in the nuclear jet. 

- The last two mode changes are seen to occur shortly before 
major radio outbursts. 

- Remarkably, no new radio knot/component was observed 
during the entire time span (1995.27 to 2008.91) covered by 
the VLBI images presented in this study - the components 
in the staircase-mode are fairly unambiguously identifiable 
with the set of components witnessesd in the straight mode. 
A smooth transition was shown in this paper, which could 
be explained by a (most likely curling) jet with a changing 
viewing angle. 

- The kinematics of the VLBI knots were modeled within the 
non-ballistic motion scenario (Gong et al. 2008). 

- The long-term radio flux variability can be explained by 
Doppler boosting effects arising from the precession of the 
nuclear jet. 

- The model of non-ballistic motion of hotspots explains the 
variability of flux-density and morphology, the significant 
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morphological changes and the unusual kinematics of the 
nuclear jet in PKS0735 + 178. 
- Similar kinematic properties - fairly stationary components 
- and a correlation between periodic ridge line changes and 
the long-term radio outbursts were found by some of us for 
two other well-known blazars: 1803+784 and 0716+714. 

Acknowledgements. This research made use of data from the MOJAVE database 
that is maintained by the MOJAVE team. This research made use of data from the 
University of Michigan Radio Astronomy Observatory, which is supported by 
the National Science Foundation and by funds from the University of Michigan. 
This research made use of the NASA/IPAC Extragalactic Database (NED) which 
is operated by the Jet Propulsion Laboratory, California Institute of Technology, 
under contract with the National Aeronautics and Space Administration. The 
National Radio Astronomy Observatory is a facility of the National Science 
Foundation operated under cooperative agreement by Associated Universities, 
Inc. The research of the authors Biping Gong and J.W.Zhang is supported in part 
by the National Natural Science Foundation of China, under grands 10778712. 



References 

Abraham, R.G., McHardy, I.M., & Crawford, C.S. 1991, MNRAS, 252, 482 
Agudo, I., Gomez, J.L., Gabuzda, D.C., et al. 2006 A&A, 453, 477 
Baath, L.B. & Zhang, F.J. 1991, A&A, 243, 328 

Britzen, S., Vermeulen, R.C., Campbell, R.M., et al. 2008, A&A, 484, 119 

Britzen, S., Kudryavtseva, N.A., Witzel, A., et al. 2010, A&A, in press 

Britzen, S., Kam, V.A., Witzel, A., Agudo, I., et al. 2009, A&A, 508, 1205 

Ciprini, S., Takalo, L.O., Tosti, G, et al. 2007, A&A, 467, 465 

Cotton, W.D., Wittels, J.J., Shapiro, I.I., Marcaide, J., et al. 1980, ApJ, 238, L123 

Fan, J.H., Xie, G.Z., Lin, R., et al. 1997, A&AS, 125, 525 

Gabuzda, D.C., Wardle, J.F.C., Roberts, D.H., et al. 1994 ApJ, 435, 128 

Gomez, J.L., Guirado, J.C., Agudo, I., et al. 2001 MNRAS, 328, 873 

Gong, B.P. 2008, MNRAS, 389, 315 

Gopal-Krishna, Stalin, C.S., Sagar, R., Wiita, P.J. 2003, ApJ, 586, L25 
Gopal-Krishna, Dhurde, S., Wiita, P.J. 2004, ApJ 615, L81 
Goyal, A., Gopal-Krishna, Anupama, C, et al. 2009, MNRAS, 399, 1622 
Gupta, A.C., Fan, J.H., Bai, J.M., Wagner, S.J. 2008, AJ, 135, 1384 
Lister et al., 2009 AJ, 137, 3718 

Nakamura, M., Meier, D.L. 2009. Proceedings of: Approaching Micro- 

Arcsecond Resolution with VSOP-2: Astrophysics and Technology, Y. 

Hagiwara, E. Fomalont, M. Tsuboi, and Y. Murata (eds.), ASP Conference 

Series, Vol. 402, 357 
Piner, Bhattarai, D., Edwards, P. G, Jones, D. L. 2006, ApJ 640, 196 
Pursimo, T., Nilsson, K., Sillanpaa, A., Takalo, L.O. 1999. In: Proceedings of 

the BL Lac Phenomenon Conference, eds. L.O. Takalo & A. Sillanpaa, ASP 

Confernence Series, 159, 385 
Qian, B. & Tao, J. 2004, 2004, PASP, 1 16, 161 
Sbarufatti, B., Treves, A., Falomo, R. 2005, ApJ, 635, 173 
Scarpa, R., Urry, CM., Falomo, R., Pesce, J.E., & Treves, A. 2000, ApJ, 532, 

740 

Spergel, D.N., Verde, L., Peiris, H.V., et al. 2003, ApJS, 148, 175 
Stickel, M., Fried, J.W., & Kiihr, H. 1993, A&AS, 98, 393 



10 



S. Britzen et al.: Understanding BL Lac objects 






Fig. 6. The VLBI images superimposed with model-fitting results are shown for all epochs between 1995.27 and 2000.03. 
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Fig. 7. The VLBI images superimposed with model-fitting results are shown for all epochs between 2001.98 and 2008.91. 
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Table 3. Model-fit results 



Epoch 


Flux-density [Jy] 


core sep. [mas] 


pos. a. [deg] 


Maj. a. [mas] Id. 


1995 Apr 7 


0.956+ 0.047 


0.00+ 0.00 


0.0+ 0.0 


0.17±0.02 


k 




0.537+ 0.038 


0.63± 0.09 


80.0± 5.4 


0.37±0.04 


cl 




0.098± 0.005 


1.25± 0.10 


38.3± 2.6 


u.yo±u. iu 


Cj 




0.05 1± 0.003 


4.35± 0.06 


62.7± 2.5 


1.47±0.15 


c5 


1995 Dec 15 


0.557± 0.028 


0.00+ 0.00 


0.0± 0.0 


0.19±0.02 


k 




0.275± 0.014 


0.62± 0.05 


94.3± 2.5 


0.32±0.03 


cl 




0.196± 0.010 


0.91± 0.05 


70. 1± 2.7 


0.33±0.03 


c2 




0.099± 0.005 


1.00± 0.05 


37.7± 2.9 


n 7Q-i_n en 


Cj 




0.040+ 0.002 


2.25± 0.08 


35.1+ 3.0 


1.18+0.12 


c4 




0.018± 0.001 


4.33± 0.07 


62.6± 2.7 


0.74±0.07 


c5 


1996 Jan 19 


0.565± 0.028 


0.00+ 0.00 


0.0+ 0.0 


0.15±0.02 


k 




0.468± 0.023 


0.69+ 0.05 


83. 1± 2.6 


0.43±0.04 


cl 




0.085± 0.004 


1.12+ 0.05 


66. 1± 2.7 


n cxi 






0.093± 0.005 


1.21+ 0.06 


29.2± 2.7 


0.87±0.09 


c3 




0.042± 0.002 


4.07± 0.06 


62.3± 2.8 


1.42+0.14 


c5 


1996 Mar 22 


0.478± 0.024 


0.00± 0.00 


0.0± 0.0 


0.15±0.02 


k 




0.370± 0.019 


0.7 1± 0.05 


83.9± 2.6 


ft AA-uC\ C\A 


C 1 




0.109+ 0.005 


1.07± 0.05 


67.7± 2.7 


0.28±0.03 


c2 




0.090+ 0.006 


1.44+ 0.07 


33.8± 3.2 


0.99±0.10 


c3 




0.033± 0.002 


4.41± 0.06 


63.2+ 2.6 


1.34+0.13 


c5 


1996 May 27 


0.446± 0.022 


0.00± 0.00 


0.0± 0.0 


0.00±0.10 


k 




0.408+ 0.020 


0.77± 0.05 


81.8± 3.2 


0.38±0.04 


cl 




0.081+ 0.010 


1.24± 0.08 


54.6± 3.5 


0.48±0.05 


c2 




0.055± 0.006 


1.42+ 0.10 


23.3+ 3.5 


n 7Q_i_n cxi 


Cj 


1996 Jul 27 


0.402± 0.020 


0.00± 0.00 


0.0+ 0.0 


0.11±0.10 


k 




0.395± 0.020 


0.79+ 0.05 


80.8± 2.7 


0.48+0.05 


cl 




0.049± 0.003 


1.23+ 0.05 


55.7+ 2.8 


0.35±0.04 


c2 




0.058± 0.003 


1.34± 0.08 


20.9± 2.7 


u.oz+u.uo 


Cj 




0.019+ 0.001 


2.51+ 0.06 


45.5± 3.0 


1.12+0.11 


c4 




0.017± 0.001 


4.14+ 0.05 


61.4± 2.9 


0.52±0.05 


c5 




0.012± 0.001 


5.80± 0.08 


65.4± 2.7 


0.86±0.09 


c6 


1996 Sep 27 


0.379± 0.019 


0.00+ 0.00 


0.0+ 0.0 


0.10±0.10 


k 




0.391± 0.020 


0.83+ 0.05 


76.9± 2.4 


U. D j±V.yJD 


Cl 




0.090± 0.005 


1.38+ 0.09 


29.7± 2.5 


1.03±0.10 


c3 




0.026± 0.001 


4.24± 0.05 


63.8± 2.4 


1.07+0.11 


c5 


1996 Dec 06 


0.327+ 0.016 


0.00± 0.00 


0.0± 0.0 


0.12+0.10 


k 




0.329± 0.016 


0.84± 0.05 


75.0± 2.4 


0.62±0.06 


cl 




0.033± 0.002 


1.42+ 0.07 


17.1+ 2.5 


U. J J±U.UO 


Cj 




0.017± 0.001 


2.25± 0.04 


35.2± 2.4 


0.73±0.07 


c4 




0.021± 0.001 


4.19± 0.05 


62.7± 2.4 


0.57+0.06 


c5 


1997 Mar 10 


0.350± 0.020 


0.00± 0.00 


0.0± 0.0 


0.22±0.02 


k 




0.251+ 0.013 


0.89+ 0.13 


75.2± 5.2 


0.66±0.07 


cl 




0.053± 0.016 


1.55+ 0.10 


22.9± 6.5 


U. /D±U.Uo 


Cj 




0.01 1± 0.002 


2.29± 0.43 


50.2± 8.1 


0.39±0.04 


c4 




0.033± 0.006 


3.96± 0.23 


64.5± 2.7 


1.43±0.14 


c5 


1998 Oct 30 


0.503± 0.025 


0.00+ 0.00 


0.0± 0.0 


0.17±0.02 


k 




0.092± 0.005 


0.77± 0.06 


75. 1± 3.0 


0.82±0.08 


cl 




0.033± 0.002 


1.83+ 0.05 


32.4± 2.9 


1 m-i-ft 1 n 


Cj 




0.019+ 0.003 


2.89± 0.06 


54.2± 2.6 


0.90±0.09 


c4 




0.018± 0.002 


3.93+0.10 


67.3± 2.6 


1.01+0.10 


c5 


1999 Sep 1 


0.402± 0.020 


0.00± 0.00 


0.0± 0.0 


0.17±0.02 


k 




0.157± 0.008 


0.32± 0.06 


76.8± 3.6 


0.19±0.02 


ca 




0.031+ 0.006 


0.84+ 0.10 


84.7± 3.4 


0.43±0.04 


cl 




0.040± 0.002 


1.83+ 0.07 


40.3± 3.3 


1.20±0.12 


c3 




U.02/± 0.0U1 


3.t>4± U.Oo 


03.3± 3.U 


1.10+0.11 


c5 


2000 Jan 11 


0.502+ 0.025 


0.00± 0.00 


0.0+ 0.0 


0.10±0.10 


k 




0.141± 0.007 


0.4 1± 0.04 


77. 1± 2.6 


0.31±0.03 


ca 




0.045± 0.003 


0.79± 0.06 


71.7± 3.5 


0.52±0.05 


cl 




0.009± 0.002 


1.69+ 0.06 


18.0± 5.0 


0.50±0.05 


X 




0.009+ 0.001 


2.22± 0.06 


44.2± 3.1 


0.54±0.05 


c3 




0.026+ 0.001 


3.55+ 0.04 


62.1±2.6 


1.06+0.11 


c5 
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Table 4. Model-fit results continued 



Epoch Flux-density [Jy] core sep. [mas] pos. a. [deg] Maj. a. [mas] Id. 



2000 May 20 


0.558± 0.028 


0.00+ 0.00 


0.0+ 0.0 


0.14+0.01 


k 




0.147± 0.008 


0.35+ 0.05 


78.6+ 2.5 


0.31+0.03 


ca 




0.077± 0.004 


0.88+ 0.05 


76.2+ 2.5 


0.51+0.05 


cb 




0.009± 0.001 


1.72+ 0.05 


28.5+ 2.5 


0.77+0.08 


X 




O 01 R-i- O 01 O 
U.Ulol U.U1U 


9 ^zL-t- O 1 f\ 


Al O-i- A Q 


O 08-t-O 1 

u.yo±u. iu 


r'X 




0.007± 0.002 


3.16+ 0.05 


63.9+ 2.8 


0.50+0.05 


c4 




0.016± 0.001 


3.99+ 0.05 


64.2+ 2.6 


0.98+0.10 


c5 


2001 Oct 22 


0.608± 0.030 


0.00+ 0.00 


0.0+ 0.0 


0.19+0.02 


k 




0.327± 0.016 


0.41+ 0.05 


56.4+ 2.6 


0.26+0.03 


ca 




o i qs+ o 01 ^ 


U.7U+ U.UJ 


f,Q Qj- 9 7 
Do.7± Z. / 


n aa+o fi7 
u.ou±u.u / 


rh 
CD 




0.022± 0.001 


2.15+0.05 


55.3+ 2.7 


1.21+0.12 


CC 




0.018± 0.001 


3.65+ 0.05 


63.3+ 2.4 


1.09+0.11 


cd 


2001 Dec 23 


0.640± 0.032 


0.00+ 0.00 


0.0+ 0.0 


0.18+0.10 


k 




0.335± 0.017 


0.46+ 0.05 


56.0+ 2.6 


0.25+0.03 


ca 




o 9s i -+- n m t, 


n Qi + n os 


fn 1+ 9 q 

U / . / IE Z.7 


78+0 08 


cu 




0.015± 0.001 


2.33+ 0.06 


55.0+ 2.8 


0.98+0.10 


CC 




0.017± 0.001 


3.52+ 0.05 


65.9+ 2.6 


1.26+0.13 


cd 


2002 Nov 23 


0.702± 0.035 


0.00+ 0.00 


0.0+ 0.0 


0.09+0.10 


k 




O 11 9-t- O 01 A 


o n_i- o os 
u.jj± yj.yjj 


AS O-t- 9 f\ 
OO.Ux z.o 


O 1 Q-t-O 09 


ca. 




0.108± 0.005 


0.84+ 0.05 


53.2+ 2.7 


0.30+0.03 


cb 




0.214± 0.011 


1.21+0.05 


63.9+ 2.6 


0.68+0.07 


X 




0.065± 0.004 


2.23+ 0.08 


62.6+ 3.4 


1.63+0.16 


CC 


2003 Aug 28 


0.640+ 0.032 


0.00+ 0.00 


0.0+ 0.0 


0.16+0.02 


k 




O 1 78-*- O OOQ 

u. i / o± u.uuy 


O S7-i- O OS 
U.jZ± U.Uj 


Aa- Q 


O ^1 4-0 0^ 


ca. 




0.224± 0.030 


0.94+ 0.06 


55.8+ 2.9 


0.56+0.06 


cb 




0.169± 0.027 


2.20+ 0.07 


60.9+ 2.9 


1.61+0.16 


CC 


2004 Aug 9 


0.523± 0.026 


0.00+ 0.00 


0.0+ 0.0 


0.17+0.02 


k 




0.103± 0.008 


0.52+ 0.06 


62.8+ 4.7 


0.22+0.02 


ca 




O 1 O OOQ 

u.io i f± u.uuy 


1 97-t- o 02 
i .z / ± u.uo 


Al S-t- zL 1 


O 71 -t-0 07 


rh 
CD 




0.101± 0.008 


2.52+ 0.09 


52.4+ 3.2 


1.54+0.15 


CC 




0.060± 0.003 


3.33+ 0.06 


65.9+ 3.0 


0.95+0.10 


cd 


2004 Dec 2 


0.522± 0.026 


0.00+ 0.00 


0.0+ 0.0 


0.17+0.02 


k 




0.140± 0.011 


0.58+ 0.07 


62.2+ 3.4 


0.32+0.03 


ca 




n 1 1 s+ o 01 9 

U. 1 U± U.ulZ 


1 .ZO± U.UU 


69 7-1- 9 f. 

DZ. / ± Z.O 


8S4-O OQ 


CD 




0.112+ 0.014 


2.52+ 0.06 


56.6+ 4.4 


1.44+0.14 


CC 




0.057± 0.010 


3.46+ 0.09 


64.1+ 2.8 


0.97+0.10 


cd 


2006 Apr 28 


0.503± 0.025 


0.00+ 0.00 


0.0+ 0.0 


0.19+0.02 


k 




0.110+ 0.006 


0.49+ 0.05 


63.0+ 2.9 


0.36+0.04 


ca 




O 1 S^4- O OOS 
U. 1 J J± U.UUo 


1 1 04- O OS 


AO 7-i- 1 8 
OU. / ± Z.o 


O AS-i-0 07 
U.OJ±U.U / 


rh 
CD 




0.026+ 0.001 


2.25+ 0.05 


51.7+ 2.8 


1.58+0.16 


CC 




0.105+ 0.005 


3.41+ 0.05 


62.7+ 2.8 


1.30+0.13 


cd 


2007 Mar 2 


0.408+ 0.020 


0.00+ 0.00 


0.0+ 0.0 


0.17+0.02 


k 




0.036+ 0.005 


0.55+ 0.06 


55.8+ 8.2 


0.04+0.10 


ca 




O 1 Al 4- O 01 f\ 
U. i L ti ± U.U 1 D 


O QS4- O OS 
\J.yJ± U.Uo 


D / .U± LJy 


O AA-i-0 07 
U.OO±U.U / 


rh 
CD 




0.055+ 0.003 


2.26+ 0.06 


53.7+ 2.9 


1.11+0.11 


CC 




0.074+ 0.004 


3.65+ 0.05 


64.2+ 2.9 


1.39+0.14 


cd 


2007 Aug 9 


0.470+ 0.024 


0.00+ 0.00 


0.0+ 0.0 


0.20+0.02 


k 




0.043+ 0.003 


0.46+ 0.06 


77.3+ 3.8 


0.30+0.03 


ca 




0.094+ 0.014 


1.06+ 0.05 


61.9+ 4.3 


0.58+0.06 


cb 




0.064+ 0.013 


2.32+ 0.13 


56.8+ 5.3 


1.48+0.15 


CC 




0.072+ 0.006 


3.54+ 0.08 


63.4+ 3.1 


1.37+0.14 


cd 


2008 Nov 26 


0.327+ 0.016 


0.00+ 0.00 


0.0+ 0.0 


0.14+0.10 


k 




0.122+ 0.006 


0.41+ 0.05 


76.0+ 3.0 


0.29+0.03 


ca 




0.072+ 0.004 


0.97+ 0.05 


67.8+ 2.9 


0.53+0.05 


cb 




0.039+ 0.002 


2.37+ 0.06 


56.2+ 2.9 


1.60+0.16 


CC 




0.075+ 0.004 


3.65+ 0.05 


63.7+ 2.9 


1.42+0.14 


cd 



